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Abstract 

We discuss semileptonic B decays of the form B J/'f ei^ X as possible probes 
of intrinsic charm. We calculate the leading order perturbative contribution to the 
process B ^ J/^ e~Ve X and find it to be unobservably small, with a branching 
ratio ~ 10-1°. We propose a modified spectator model to estimate the intrinsic charm 
contribution and find that it can be significantly larger, with a branching ratio for 
B (cc) e~lj(, X as large as 5 x 10^^. We show that the process could be observed at 
these levels by CDF assuming a Run II integrated luminosity of 15 fb-^, making this 
a useful reaction to probe the idea of intrinsic charm. 



1 Introduction 



Interesting possibilities open up if we consider wavefunctions of hadronic bound states that 
contain light-cone Fock states of arbitrarily high particle number. The higher light-cone 
Fock components arise as quantum fluctuations suppressed by M^, where M is the mass of 
the fluctuation. The cc component of hadronic Fock states is referred to as intrinsic charm 
(IC) i- 

The possibility of detecting signatures of IC in 5-meson physics has been recently 
discussed in Ref. [H. The authors propose IC as the explanation for a "slow" (low momen- 
tum) bump in the inclusive B J/"^ X spectrum and the softness of the J/\E' spectrum in 
T — s> J/\l' X decays, contrary to expectations from the color octet mechanism. IC appearing 
in higher light-cone Fock components of the i?-meson wavefunction can manifest itself in two 
possible ways. It can operate virtually, in a mediation role, and affect decay processes by 
providing additional channels for the weak interactions. In this way one may enhance some 
CKM suppressed B decays. This has been discussed in Ref. using a phenomenological 
mixing angles approach, and more recently in Ref. 0] . Another instance of this mediation is 
given by the hypothesis that the IC component of p mesons is able to explain the "pvr puz- 
zle" . IC can also manifest itself in processes in which the (intrinsic) cc produces charmed 
hadrons in the final state. An example of this is discussed in Ref. Q to account for charm 
production in deep inelastic scattering. Here we consider another example of the latter by 
looking for J/\E' production from a cc component in the B meson wavefunction. 

To determine the process that is best suited for our purpose, we first consider purely 
hadronic decays of the B meson with just one cc pair in the final state. It is then easy to see 
from Fig. [|, that hadronic channels with a cc Xu,d,s quark content are also produced at tree 
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level via V^^, with no need for IC. We conclude that these processes are ill-suited to study 
the effects of IC, which would be swamped by a formidable background. Next, we consider 
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Figure 1: Perturbative contribution to the decay — > J/^! X. (a) is a standard tree-level 
graph, (b) is an annihilation graph involving IC in the initial state. An IC exchange graph 
is also allowed. 

a process with three charm quarks (or antiquarks) in the final state such as i? — J/"^D^*^ 
or B J/'^D^*^X. The former has been discussed in Ref. and in the latter the phase 
space available closes rapidly, suppressing the rates to unobservable levels. 

This leads us to concentrate in semileptonic modes with a J/\& in the final state. These 
modes have the additional advantage of presenting an extremely clean experimental signature 
containing three charged leptons. The simplest mode of this type is i? J/\E' e~z7eX.[] 
Detectors such as CDF operating at hadron machines, or experiments at B factories such 
as BABAR at SLAC should be able to identify the three outgoing charged leptons easily. In 
fact, CDF will be looking for the mode Be ^ J/"^ e~Ve which has the same topology (and 

which would constitute a background to our process). It would therefore be natural for them 

^For definiteness we concentrate here on channels with a J/'^ in the final state, but our analysis can also 
be applied to channels with other cc resonances or with DD pairs instead. 
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to search for i? — J/\l/ e Ue ^ simultaneously. 

This paper is organized as follows: in Section ^ we calculate the perturbative QCD 
contributions to semileptonic processes such as J/\E' eu X and show that they are 

unobservably small. In Section ^ we propose a modified spectator model for the computation 
of the IC contributions to these same semileptonic processes, and find that they can be 
substantially larger, with branching ratios as large as 5 x 10^^. Finally in Section |^ we 
estimate the reach of near future experiments for these reactions and show that they could 
observe the mode B~ —* J/\E' eu X if the IC component of B mesons is as large as has been 
speculated recently 0. 

2 The perturbative calculation 

In this section we compute the perturbative contribution to the decay rate of the process 
B^ ^ J/^i e~VeX from the two lowest order diagrams in Fig. |^. We present calculations 
of the decay rate in both Non-Relativistic QCD (NRQCD) |^ and in the Color Evaporation 
Model (CEM)P]. We also derive the leading order contribution to the total inclusive B-to- 
(cc) semileptonic decay rate. 

A hadronic B decay depends on the momentum scale of the meson wavefunction 
(1/fm ~ Aqcd) and is non-perturbative in principle. Nevertheless, the mass of the heavy 
quark provides a large momentum scale and to the extent that rrij^^ ~ 2mc ^ Aqcd and 
— rriji^ ^ Aqcd, the cc pair production occurs at short distance and is perturbative. On 

the other hand, the formation of J/\E' mesons from the initially compact cc pairs takes a long 
time and is non-perturbative. The quantum interference between the production of the cc 
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Figure 2: Graphs describing the perturbative contribution to the decay B — J/\E' e z/X. 

pair and the formation of J/"^ mesons should be suppressed by powers of the ratios of the 
typical momentum scales of these two processes, such as ^qcd/^'j/ip ^qcb/ i^b~^J/ipy ■ 
Therefore, if we neglect these power corrections, the decay rate can be factorized as 



X 5(f + ih - k,f) F[,,]„^j/^{q^), (1) 

where J2n ^^ns over all possible singlet and octet color states, and = A;^ + is the total 

momentum of the produced cc pair. In Eq. (|l|), = (2A;c)^ is the square of the relative 
momentum between the c and the c in their rest frame. Treating the c or c approximately 
as being on their mass-shell, = = m^, = {kc — fcg)^ = ~ 4:iTi^- The non- 
perturbative quantity F[cc]„^j/^(g^) denotes the transition probability for a [cc] state with 
relative momentum squared to evolve into a physical J ftp meson. Different choices of 
-^[cc]„^j/V'(9'^) Isad to different models of J/"^ formation. 
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2.1 Nonrelativistic QCD 



The NRQCD formalism for quarkonium production exploits the fact that the heavy quark ve- 
locity Vc is much less than one, or equivalently that the heavy quark three-momentum is much 

smaller than its mass {\kc\ -C rric) [Q- By expanding the g^-dependence of dV -,[cc]nX I d^Q 
around = in Eq. (|l]), we effectively derive the leading term in the NRQCD formalism: 

T{B- ^ J/m e-u,X) = Y.^{B- ^ [cc]ne-V,X) {Oi'^) . (2) 

n 

The constants {O:^^) oc / dq^Fycg^r.^j/ipio^) are matrix elements of local operators in coor- 
dinate space, often called NRQCD matrix elements. They contain all the non-perturbative 
information on the J/\E' formation from the cc pair. According to the NRQCD formalism 
for heavy quarkonium production, the color of the gluon in Fig. |^ is neutralized by non- 
perturbative soft gluons in the production process, which are not explicitly shown. For the 
lowest order diagrams in Fig. ^ only a spin-1 and color-octet cc state is able to produce a 

J/\&, and the corresponding matrix element is {O'l^^ i^Si)) 0. 

To the extent that the J/\E' mass is large, the partonic process in the decay amplitudes 
of \a^^e~Vf,X in Fig. ^ takes place at short distance ~ More precisely, the 

gluon is off-shell by the J/\l' mass, and the denominator of the t- or u-channel fermion 
propagator is at least 5 GeV^. Therefore, the decay in Fig. ^ can be factorized into a 
hadronic matrix element for B — > fo-u, and a partonic process hu \cc\^e~Vf,- In order to 
achieve this factorization we formally write the decay amplitudes of -B ^ g*W* in Fig. Q as: 

Mr = / ^ Tr [%ipB, p,) . nrip,, P^, q)], ^ = 1, 2 , (3) 

where Tg and H represent the matrix elements oi B ^ hu and the partonic process hu — * 



g*W* respectively. The partonic momenta satisfy — Pb ^ Phi ^'^^ leptonic 
momenta is g = + g2 = Q'e + 9f- 

By decoupling the spinor trace between the matrix element and the partonic process, 
and boosting to an infinite momentum frame, the decay amplitudes in Eq. (|^) can be 
approximated as (Mf^ and refer to the contributions from the direct and crossed graphs, 
respectively) , 

Mr ^ jdx T^{x) HZip, = xp^, p^, q) , (4) 
where the matrix element is defined by 

T^ix) = J exp{^xpt y') { | 11(0)^1, h{y-) \ 5") , (5) 

and the partonic process H is computed extracting the axial- vector component of the lowest 
order amplitudes hu g*W* . 

Eq. is given in terms of light-cone coordinates: the ± components of any four- 
vector k^^ are defined as = {k^ ± k^)/y/2. We also write = (j)^, m|/(2pj), 0±j using 
the boost-invariant light-cone wavefunction language. 

Due to the heavy quark mass, the typical virtuality of b and u in the B meson should 
be much smaller than the large momentum scale in H, and therefore we can approximate 
the parton momenta and in as 

Pb = xPB , Pu= i^-x) Pb , with x = — . (6) 

This approximation corresponds to dropping power corrections suppressed by the hard mo- 
mentum scale of the partonic process. Fixing x in H°'J^ we decouple the x integration in 
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Eq. m to find 



dx T^{x) = ( I u{0)r% m I B-) = if J, p% 



(7) 



For numerical estimates we use = 0.177 GeV [10, 11||. Substituting Eq. into Eq. 
we obtain a factorized expression for the decay amplitude of B ^ g*W*, 



Mr ~ {ifs Pb) H^J" (pfe = xPb, Pii,, q 



The amplitudes in Eq. (H) can thus be written as 



Ml 



9s9w 



B 



2^/2 {Pm-Puf - 



-V r 

2 ub 



e^'^^PB^ip^ - P^) 



+ Pb^P^ - PuT + Pb{p^ - PuY - Pb ■ (P* - Pu) 9 



(9) 



and (for the crossed diagram), 



A/f'^A' _ f 9s9w 1 



2^2 {q-Pu) 



e^'^-^p^Ap^ - q), 



+ pb{-Pu + qT + + qf - pb ■ i-Pu + q)9 



(10) 



We denote the strong and weak couphng constants by and respectively. 

The W* coupling to the lepton current is evaluated in the usual way. After squaring 
this contributes a factor 



L"^ — gw I ~^ 2~ 



'111 



^Ib = 0.177 GeV corresponds to the usual annihilation of a bu pair in the color singlet configuration. 
In our case the bu pair is in a color octet configuration and therefore we expect that our fg is actually 
smaller than 0.177. Our results for the perturbative production of J/*!/ in semileptonic B decays can thus 
be considered to be an upper bound. 
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to the decay rate. 

Our remaining task is to parametrize the production of a J/\E'(^S'i) state from the cc 
pair. We do this with projector techniques and compute the overall color factor below (in 
Eq. (|15D ID). We write 

JM) = E 1 ^A{-Ws%) n?] ef{p^) ■ (\\ Tr[fe7.) nf] e<'\p^) , (12) 
A VP*/ VP*/ 

with the projector defined as 

ni = {^^ - m^) 7^ (i** + m^) . (13) 



The total contribution to the decay rate from the cc line is then 



1 \' 1 



\pi J ^ ^ ^ J 



and the overall color factor is 



E nT^Ci] Tt[T^C^] = -, Ci = V2T\ (15) 

c A.BC 

where T^, A = 1...8, are the generators of the color group. 

As shown in Eq. (H), the rate of the process J e'V^X includes the NRQCD 

matrix element ((9g^*(^5'i)) defined in |12|, [1^, which represents the probability for an 



almost point-like cc pair in a color-octet '^5'i state to bind and form a J/\I^. The matrix 
element is nonperturbative and has to be extracted from experimental data, and its numerical 
value is I, |12|: 

{Oi'^'i^Si)) = 1.19 X 10-2 GeV^ . (16) 



As a final step to obtain the decay rate we need to divide by the number of color and 
polarization states of the cc{^Si) produced, N^oiNpoi = 24, as explained in Ref. 



Collecting our results, we obtain for the matrix element squared 



\M\' = (Mr + [Mt + Mr ) JM) ^ 



rl3u 



1 {oi^^es,)) 



2 N,oiN, 



pol 



(17) 



which includes the color factor 1/2 obtained in Eq. (|T^). In our numerical results we use the 
boson masses ITl] = 5.279 GeV, = 3.097 GeV, and the current quark mass = 
4.9 GeV, together with \V^f,\ = 0.0035 |T§, as{2m^) = 0.266 [|T|] and = 1.653 x lO^^^ 



sec WM. We find a branching ratio: 



BR{B- J/m e'VeX) = 9 X 10 
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'1^ 



The process in Fig. ^ which we have calculated is the simplest one giving rise to 
B J/\l' e^Vf. X. However, it suffers from a suppression factor associated with the 
hu annihilation. For this reason it is possible that a process involving more quarks but no 
annihilation might actually be larger. To investigate this possibility we now estimate the 
rate of the process B ^ J /'^ e^z/g -^m,c as in Fig. |^. 

We adopt once more the NRQCD factorization form of Eq. (|^) to calculate separately 
the matrix elements of 6 — > g g^e'Ty^. [q = u,c), and for the conversion of the cc pair into a 
J/^. We denote by iff and the amplitudes of the partonic processes where the gluon 
is emitted from the b and q lines respectively. The squared partonic matrix element is then 



M"^'^'' = Tr {H^'' + i/7)(fff" + 



(19) 
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Figure 3: Spectator graph analogous to Fig. ^ for i? ^ J/\E' e v X^^^- The gluon propagator 
may also be attached to the outgoing quark u or c. 



with 



Hr = (p^ - pI)2 _ ^2 - 75)( A- Ml'uM 



Hr = 'fjl _ _ KA(P.)7^( A- ^)7"(1 - 75)^.(P.) • (20) 

The leptonic and hadronic cc parts and the NRQCD matrix element are the same as for 
B — > e^UeX, and have been given in Eqs. {^^, (p!4D, and (|16[), respectively. The 
overall color factor is different, 

-1 E Tr[r^C8^] Tt[T^C^] E (t^) . . (t^) . . = I , = . (21) 

^^c ABC ij 

All these factors are combined into the squared matrix element, 



(22) 
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Using the masses already defined plus = 1.25 GeV and \Vf^^\ = 0.041, and a b quark 



lifetime = 1.6 x 10 sec we find: 



These results show that the diagrams in Fig. |^ generate the leading perturbative contribution 
to the process B ^ J/"^ e^Ve Xq, which is unobservably small. 

2.2 Total inclusive B-to-cc semileptonic decay 

The inclusive semileptonic rate for cc e^Vf, X is obtained from the diagrams in Fig. ^ 

with some simple replacements. The calculation proceeds as in Section |2]1| up until Eq. (]TT|). 



After that we simply couple the g* to a cc line in the usual way, obtaining 

where Pi and P2 are the momenta of c and c, respectively, and p = Pi + P2- The overall color 
factor now takes into account both the color octet and the color singlet configurations of 
the incoming quarks b and u, and can be calculated to be 2/9. The squared matrix element 
becomes 



\M\' = (Mr + + Mf L^, J',M,P2) 7T > (25) 



2 
9 

where and L^^ are the same quantities already used in Eq. ([T7|) . 

We use the masses and CKM parameters defined previously, but we require the in- 
variant mass of the cc pair to be larger than 3 GeV in order to produce a physical cc meson. 
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We obtain a branching ratio 

BR{B- cc e'V^ X) = 4.2 x 10"^°. (26) 

The differential branching ratio is plotted in Fig. ^ as a function of the cc pair invariant mass. 
From Fig. ^ we can see that about 90% of the time the cc pair will become a charmonium 



dBR(B"^cc e~i/)/dm^- (GeV~^) 




3.0 3.5 4.0 4.5 5.0 

m,^ (GeV) 



Figure 4: The differential branching ratio of i? —> cc e z/ X is plotted vs. the invariant mass 
rricc of the cc pair. The vertical dashed line shows the energy threshold for DD production. 

resonance and only about 10% of the time it will lead to DD production. 
2.3 Color evaporation model 

According to the Color Evaporation model of J/\E' production, all the cc pairs with invariant 
mass below open-charm threshold should become charmonia by radiating soft gluons to 
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neutralize the color and other quantum numbers The formation process from the cc 
pairs to J/\l/ mesons is assumed to be independent of the color and invariant mass of the 
pairs. This is equivalent to choosing 



(27) 



in Eq. (|T]). Here 4m|, represents the open-charm threshold, and /j^^ is a universal constant 

encoding the relative probability for a cc pair to become a J/\l/ instead of other charmonium 
states, and needs to be fixed by data. Substituting Eq. ( P7D into Eq. (|^), we derive the 
formula for the decay rate in the Color Evaporation model. 



B--^J/il) euX 



J 2mc 



drrirr. 



(28) 



where rricc is the invariant mass of the cc pair. Using the results of our previous section 
shown in Fig. ^, this becomes 



BR{B- J/ijj evX) = fji^ 3.8 x 10 
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(29) 



Comparing Eq. (|29|) with Eq. ([T8|) we see that the CEM model result would agree with the 



NRQCD result if f jj^ 
in Ref. HI. 



1/4.3. This fits very well with the expectation that 1// 
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Our results in Eqs. (|T8D and (^31) show that the perturbative QCD contributions to 
the processes J/\E' e~VeX and B J/^ e~Ve Xu,c are too small to be observed by 

near-future experiments. Qualitatively, these results are so small because they involve two 
off-shell propagators as well as suppressions from powers of and from the octet matrix 

element {Og^^ {^Si)). This leaves open the possibility of significant enhancements due to IC 
which we explore in the following section. 
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3 The intrinsic charm contribution 



The mechanism by which IC in a -B meson leads to the decay J/'^ e~T7f.X is dif- 

ferent from the perturbative case. The two (anti)charm quarks are now part of the initial 
wavefunction as shown in Fig. ^. In this picture a fluctuation in an ordinary B~ meson 

b 



B 




Figure 5: IC contribution to the decay B ^ J e u X. 

produces a nearly on-shell cc pair which is part of the initial wavefunction. There is a small 
probability for this fluctuation to occur |^ , but once it happens, it is possible for the off-shell 
hu quarks to annihilate leaving a J/\I' behind. In this section we attempt to estimate the 
rate of this process in the following way. We start from a B~ meson wavefunction in the 
higher Fock state \bucc) with a certain momentum distribution as proposed in the IC picture 
of Brodsky et al. We then estimate the annihilation rate of the hu components treating 
the cc as spectators. This annihilation is negligible unless the hu system is in a J = 1 state, 
and in this case we think of it as an off-shell B* . As we will see, the invariant mass of the 
remaining cc pair is dominantly in the charmonium resonance region, so we assume that a 
certain fraction produces a J/\I' as in the CEM model discussed in the previous section. 
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The momentum fraction distribution of the \bucc) hght-cone Fock component of B~ 
imphes that IC carries large momentum fractions, unhke the small x of extrinsic quarks from 
gluon splitting. It is given by [jl[: 

dPFc ^ S{1-EUx,) 



dxi ■ ■ ■ dxi {m\ — m\lxu — ^l/^b ~ ^l/^c ~ ^11 ^c)"^ 

In Eq. (^) the x^ are the light-cone momentum fractions + kr')i/{p^Q and rhi 

= (mf + {k\)iY/'^ are the transverse masses of each constituent as in Refs. DH. In our 
calculation we use the numerical values mb = 5 GeV, niu^ = 0.45 GeV and fric = 1.8 GeV 
[jl], [l^. Integrating Eq. (^) over all possible momenum fractions gives the probability Rjc of 



having the B meson in an IC configuration. In terms of this Rjc we obtain the four-particle 
normalization factor A^4 = Rjc x 10^ GeV^. It has been suggested in Ref. [Q that IC effects 
in B mesons may be four times larger than in the proton, where it is estimated that they 
can occur at the 1% level @. In particular, this means that Rjc could be as large as 0.04 
according to Ref. [Q. 

We are interested in a process in which the cc is in the particular configuration ^Si 
color singlet state. There is, of course, a model dependent combinatorial factor associated 
with this. Because we lack a detailed dynamical model for IC, we will instead resort to a 
phase-space model similar in spirit to the GEM discussed in the previous section. 

We begin our estimate of the amplitude by treating the charm (anti) quarks as spec- 
tators. The wavefunction is then the product of a (cc) state that generates the J/\E' and 
a virtual (bu) state with the quantum numbers of a B*^ (other configurations are helicity 
suppressed) which decays into leptons. This virtual {bu) pair possesses an invariant mass 
squared = <f = (a^b + XuY'''^'b- The amplitude squared reduces to the product of a lep- 
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tonic part L^^ which is the same as that in Eq. (|TT]) and the hadronic b-u hne contribution. 
Combining them we find: 

\M\' = -/l^^l^f^lKJV^. = ^fU'Gl |K.r • (31) 

It is evident from this result that we treat the bu annihilation as the leptonic decay of an 
off-shell B*. For our numerical estimates we take = f^. 



It remains to integrate Eq. (|3TD over the phase space available to the lepton pair 
(determined by m^*) and to convolute this result with the momentum distribution in the 



B meson wavefunction, Eq. (|30D . In this convolution we impose kinematic constraints 
consistent with "liberating" the IC pair which becomes a charmonium resonance or a DD 
pair. 

Numerically, all this results in the differential decay rate shown in Fig. ^, where we 
have used Rjc = 0.04 for illustration. The integrated rate leads to the inclusive branching 
ratio, 

BRjciB- ^ (cc) e-VeX) = 5 x 10"^ x ^ . (32) 
If instead we only integrate over the charmonium resonance region, we obtain the estimate, 

BRjciB- ^ J/^ e-V, X) = 4 X 10"^ x ^ . (33) 

Only a fraction of this result corresponds to J/"^, and if we use our perturbative calcula- 
tions as guidance, this fraction could be of order /j/^ ~ 1/4.3. However, as a signal for IC 
any charmonium resonance that can be reconstructed (or for that matter any DD channel) 
will serve. Integrating the distribution in Fig. ^ above open-charm threshold leads to the 
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Figure 6: The IC contribution to the differential branching ratio of -B cc e z/ X is plotted 
vs. the invariant mass m^c of the cc pair. The vertical dashed line shows the energy threshold 
for DD production. 

estimate BRj^(5~ DD e'VeX) = 6 x 10"^ Ric/^M. If the fraction of IC in the B 
meson can really be as large as Ric = 0.04, we see that semileptonic B decay constitutes a 
powerful probe of these ideas. Our result, Eq. (|33|) , is three orders of magnitude larger than 
conventional, perturbative, mechanisms for producing charm pairs in semileptonic B meson 
decay. This large enhancement can be understood qualitatively because the PQCD mecha- 
nism is suppressed with respect to the IC mechanism by two powers of ag, two propagators 
with an m^, mass scale, and the matrix element for J/\E' formation from a color octet cc. 
At the same time the only relative suppression in the IC mechanism is the probability of 
finding the B wave-function in the IC configuration, Ric- Roughly, the ratio of PQCD to 
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IC induced rates goes as al{{Oi'^ Si)) /ml) / Ric ~ 10-=^. 

We can perform a similar calculation for the decay B ^ J /'^ e^v^, Xu^c-, depicted in 
Fig. 0. This process has at least one more particle in the final state and we might expect 
it to be suppressed by phase space with respect to i? ^ J/\& e^V^X. However, it does not 
suffer from the suppression factor associated with the weak annihilation of the hu pair and for 
this reason we estimate it here. Our treatment of IC is similar to the previous channel. We 
calculate the semileptonic decay of an off-shell 6-quark into three particles, h — e^Vf. {u or c) 
and treat the remainder three quarks in the B wavefunction as spectators. 

As in our previous example, we picture a nearly on-shell cc pair arising as a fluctuation 
in the S-meson wavefunction. The semileptonic decay of the 6-quark is standard, the only 
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Figure 7: IC contribution to the three-body spectator semileptonic decays h ^ J e v^, 
difference here being that the 6-quark is off-shell and carries a fraction x\, of the B meson 
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momentum. The partonic decay rate is thus (in terms of ml = x^m^): 

, GllKfrnf (m\ 
r(6-.e-..)«^i|L_i/y . q = u,c, (34) 

where for the case of a V^^ transition we also include the usual kinematic factor 

This factor is approximately equal to 1 for the V^^ transition. This result is now folded with 
the momentum distribution in the B wavefunction, Eq. (30), with kinematic constraints to 



once again "liberate" the IC pair. In spite of the CKM suppression Vub/Vcb, the Vub mode can 
be as large or larger than the Vd, mode because there is a very large phase-space suppression 
producing three charm quarks. This suppression is so strong that the result is very sensitive 
to the value of mc. Using = 1.25 GeV and a b quark lifetime = 1.6 xlO"^^ sec, we 
estimate 

BR,^(5 J/vf eu, ^ I i ^ Zl ^ ^^^/O-O^ f ^ ^ = ^ ' (36) 
' 'i' |4 X 10"^ X /?/c/0.04 forg = M. ^ ' 

If instead we use = 1.6 GeV we obtain BRjc{b J/^ e'Ve X J = 4 x 10"^ x Rjc/OM. 
From this we conclude that the V^b mode is a small correction to our previous process so 
that — s> J/\E' e^i? is dominated by the annihilation mode B~ J/^ e^u X. We also 
conclude that the channel with three charm quarks in the final state is smaller in spite of 
the CKM advantage it enjoys. 



4 Discussion and conclusions 



Up to now we have discussed the production of a cc state like J/\E' in a semileptonic B 
decay. We have shown that the perturbative QCD contributions to this process are extremely 
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small, giving rise to branching ratios at the 10 level. Similar conclusions apply to the 
perturbative production of other charmonium resonances or DD pairs. We have also shown 
that there can be significantly larger contributions from IC components of the 5-meson 
wavef unction. 

It is simple to understand why the IC contribution to our processes is much larger than 
the perturbative contribution. The perturbative contribution suffers from several suppression 
factors associated with the exchange of a perturbative gluon and the color rearrangements 
that are needed to convert two color octet qq pairs into color singlet mesons. The IC 
contribution bypasses all this because all necessary color exchanges and rearrangements take 
place in a nonperturbative way inside the 5-meson wavefunction. The only suppression for 
this mechanism is then the small probability of finding the B meson in the IC configuration. 

Adopting a simple spectator model and accepting recent suggestions in the literature 
as to the possible size of the IC component, we have estimated that (cc) e^vX could 

have a branching ratio as large as 5 x 10~^. We now show that at this level, the process 
B^ (cc) e^vX could be observed in the near future. 

We start by estimating the sensitivity of CDF to the process S — J/\E' evX in Run 
II by comparing it to B^ —>■ J/\I' lu, which they expect to see and which has the same final 
state. Assuming 2 fb~^ and 15 fb~^ respectively for the integrated luminosities during Run 
Ila and Run Ila+Run lib, and using |]18| 



BR(6 ^ 



BR(5, ^ J/^ lu) 



0.13 ±0.06 



(37) 



BR(6 ^ B-) BR{B 



- J/^ K-) 



20 



leads to the estimate [| 

iV(R ^ J/^ lu) ~ 800 = 6000 1 . (3^ 

^ ^ ^ 2fb"i 15fb~^ ^ 

From this estimate, with / = /i, e, the number of events for — J/\E' /z/X is given by 



^ ' ' 15 fb"^ 10-7 ^ ' 



Therefore, CDF could be able to use B J/\l' lu X to test IC in the upcoming Run 
II. Notice that the open-charm channels such as i? ^ D^D^ lu X could also be used. 
Our modeling of IC predicts a smaller branching ratio, of order 6 x 10^*^. However, the 
identification of charged D mesons could be up to 10 times more efficient than that of J/"^ 
since the latter is reconstructed from its leptonic decay which has a small branching ratio. 



S(J/^ ^ + e+e-) ^ 0.12 

A similar projection for BABAR can be based on the decay B^ D* lu. Given the 
7517 ± 104 reconstructed events for B^ D* lu with a total integrated luminosity of 9.7 



fb ^ [T^, and a projected total luminosity of ~ 542 fb ^ for year 2006, we have: 



f r 

N(B° lu) ~ 4.2 X 10^ — r , (40) 

^ ' 542 fb"^ ' ^ ^ 



and from this we can estimate 



N{B J/m euX) ~ BR(J/* ^ ee + /x/i) x A^(5" ^ D* lu) x 



, BRjc{B^J/^euX) 



BR(50 ^ D* lu) 



10-6 ^ ^ 



^We are grateful to T. LeCompte and M. Tanaka for providing us with this estimate. 
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We have included a factor BR(J/\1/ — ee + fxfi) to roughly account for the lower re- 
constructing efficiency of J/^ compared to D*. From these numbers we conclude that 
BRjciB -> J/^ euX) could be probed at the 10"'^ level at CDF and at the 10"^ level at 
BABAR. Similar considerations can be applied to Belle. 

Near future experiments could use semileptonic S-meson decays to probe the ideas 
of IC at an interesting level. Their study would complement the suggestion of Ref. and 
can place limits on the IC content of B mesons. 
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